Low- Temperature Specific Heat of an Extreme Type-II Superconductor at High 

Magnetic Fields 
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We present a detailed study of the quasiparticle contribution to the low-temperature specific heat 
of an extreme type-II superconductor at high magnetic fields. Within a T-matrix approximation for 
the self-energies in the mixed state of a homogeneous superconductor, the electronic specific heat 
is a linear function of temperature with a linear- T coefficient Js(H) being a nonlinear function of 
magnetic field H . In the range of magnetic fields H > (0.15 — 0.2)H C 2 where our theory is applicable, 
the calculated js(H) closely resembles the experimental data for the borocarbide superconductor 
YNi 2 B 2 C. 
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In the last two decades and particularly since the 
discovery of high temperature superconductors (HTS) 
almost all of the superconducting systems that hold 
the greatest promise for practical application are of the 
extreme-type-II variety. These materials are character- 
ized by high transition temperatures (T c ), high upper 
critical fields (H C 2), and can be defined as materials in 
which the semiclassical H C 2(0) in units of Tesla becomes 
comparable to, or even larger than, T c in units of Kelvin. 
In such systems, at low temperature and high magnetic 
fields near and around the semiclassical H c2 (0), Landau 
level quantization of electronic energies within the super- 
conducting state is well defined and has to be included 
in the description of the superconducting instability 1,2 . 
Such a regime in which the Landau level structure is 
well resolved , i.e. cyclotron energy huj c > A(T, H),T, T 
(where A(T, H) is the BCS gap, and V is the scatter- 
ing rate due to disorder), represents a large portion of 
the H — T phase diagram of an intrinsically extreme 
type-II superconductor. This region can extend down 
to magnetic fields as low as ~ (0.2 — 0.5)H c2 (0) and tem- 
peratures as high as ~ 0.3T c o. In contrast, the size of 
the similar region in an ordinary type-II superconduc- 
tor (such as Nb) is expected to be negligible and con- 
fined to an immediate vicinity of H C 2(0)£ Inside this 
high-field, low-temperature regime the superconducting 
state fundamentally differs from the familiar low-field 
mixed phase of the Abrikosov-Gorkov theory<4, primar- 
ily by the appearance of gapless quasiparticle excitations 
at the Fermi surface. These gapless excitations reflect a 
coherent quasiparticle propagation over many unit cells 
of a closely packed vortex lattice with fully overlapping 
vortex cores^. 

The presence of such low-lying quasiparticle excita- 
tions makes an s-wave, "conventional" superconductor in 
a high magnetic field somewhat similar to an anisotropic, 
"unconventional" superconductor with nodes in the gap. 
In the low-temperature and high-field regime, however, 
the nodes in the gap reflect the center- of-mass motion of 



the Cooper pairs in the magnetic field, in contrast to d- 
wave superconducting cuprates where such nodes are due 
to the relative orbital motion. This gapless behavior in 
three-dimensional systems is not restricted to fields very 
close to H C 2 but rather persists to a surprisingly low mag- 
netic fields H* ~ (0.2-0.5)i?c2££<£ Below H* gaps start 
opening up in the quasiparticle spectrum and the sys- 
tem eventually reaches the low-field regime of localized 
states in the cores of isolated, well-separated vortices 8 . 
Recently, an extensive numerical calculation of quasipar- 
ticle excitations in the mixed state for both s-wave and 
d-wave superconductors was performed 9 and it was found 
that for fields H > 0.5H C 2 no qualitative difference in be- 
havior can be seen between s- and <i-wave cases. They are 
both characterized by coherent low-lying Landau level- 
like quasiparticles excitations. However, a marked differ- 
ence appears at lower fields H <C 0.5H c2 , where an s- 
wave superconductor is clearly in the regime of localized, 
bound vortex core states while a d-wave system still ex- 
hibits the extended nature of low- lying quasiparticle exci- 
tations as predicted by Franz and TesanovioiS. Further- 
more, it was shown that the high-field gapless character 
of the excitation spectrum is not destroyed by a moder- 
ate level of nonmagnetic impurities present in either dirty 
homogeneous superconductor— or dirty inhomogeneous 
superconducting systems^. 

The strongest evidence for Landau level quanti- 
zation within the superconducting state comes from 
the experimental observation of de Haas-van Alphcn 
(dHvA) oscillations in various A-15 and borocarbide 
superconductors 3 :. The persistence of the dHvA signal 
deep within the mixed state of these three-dimensional 
extreme type-II systems can be attributed to the pres- 
ence of a small portion of the Fermi surface containing 
gapless quasiparticle excitations, surrounded by regions 
where the gap is largei&i^ii At the same time, care- 
ful measurements of thermal properties (i.e. thermal 
transport and/or specific heat) at low temperatures and 
high magnetic fields are also expected to reveal the novel 
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gapless behavior in these systems. The presence of ex- 
tended gapless quasiparticle states at low temperatures 
should lead to qualitatively different thermal behavior 
than those found in an s-wave superconductor at low 
fields, where the number of quasiparticles excited above 
the gap is exponentially small and the only contribution 
might come from the bound states localized in the vor- 
tex cores. Recently, we studied the quasiparticle con- 
tribution to the thermal conductivities Kij(Tt,T) of an 
extreme type-II superconductor placed in magnetic field 
H such that H c \ « H < H c2 . We examined the trans- 
port coefficients Kij /T in the limit of f2 — ► and T — ► 
and found that there was considerable enhancement of 
thermal transport in the mixed state of an s-wave super- 
conductor due to the creation of gapless excitations in the 
magnetic field. This is in contrast to the zero field ther- 
mal transport which is exponentially small for an s-wave 
superconductor with no nodes in the gap— The agree- 
ment of our theoretical curves with the experimental data 
for the borocarbide superconductor LUN12B2C and A-15 
superconductor V3S1 by Boaknin et a/ i 17 i 18 is very good 
over a wide range of fields used in the experiments. 

The low-temperature electronic specific heat C(T,H) 
is yet another probe of the quasiparticle excitations in 
the mixed state of a superconductor. In a fully gapped 
s-wave superconductor at low magnetic fields there is 
an exponentially small contribution to C(T, H) at low 
temperatures and the only significant contribution to 
C(T, H) comes from the quasiparticles localized near the 
vortex axisiS. Assuming that the vortex core can be 
approximated as a "normal" metal embedded in a su- 
perconducting medium, this contribution is then propor- 
tional to the quasiparticle density of states which is finite 
and approximately equal to its normal state value. From 
here it follows that C(T, H) varies linearly with T and 
j s (H) = C{T,H)/T is proportional to H as T -> 0. 
On the other hand, it was predicted that in unconvential 
d-wave superconductors the density of states and there- 
fore the linear-T specific heat coefficient varies as \[H 
at low fields H > H&. This field dependence is a con- 
sequence of the delocalized quasiparticles that can move 
along the nodal directions of the order parameter—. Ex- 
periments on HTS materials YBa2Cu307_5 (Ref. 20) 
and La2_ a; Sr a ;Cu04 (Ref. 21) have confirmed this theo- 
retical prediction. The consensus has been reached that 
nonlinear field dependence of J S (H) in HTS systems is 
one of the signatures of an order parameter with d x 2_ y 2 
symmetry. However, a number of experimental studies 
do not conform to this interpretation: a nonlinear in- 
dependence of j s {H) in almost the entire regime of the 
mixed state is observed in s-wave superconductors, such 
as A-15's V 3 Si (Ref. 22) and NbSe 2 (Refs. 23 and 24) as 
well as in borocarbides superconductors LUN12B2C (Ref. 
25) and YNi 2 B 2 C (Refs. 24 and 26). It is well estab- 
lished that A-15's are fully gapped superconductors at 
zero field while there might be a significant anisotropy 
of the borocarbide's s-wave order parameter at very low 
fields H ~ iJ cl> iL2i Sonier et alM> accounted quantita- 



tively for this nonlinear behavior of J S (H) in NbSe2 by 
the expansion of the vortex cores at low fields. On bal- 
ance, it seems that the specific heat behavior of s-wave 
superconductors in the whole regime of the mixed state 
is not fully understood and merits further attention. 

The purpose of this work is to examine in detail the 
quasiparticle contribution to the low-temperature spe- 
cific heat in the mixed state of a three-dimensional s- 
wave extreme type-II superconductor starting from the 
high field limit. It was already suggested in Ref. 5 that 
in a pure superconducting system close to H C 2, the spe- 
cific heat C (H, T) should be an algebraic function of 
temperature with the power dependent on the dimen- 
sionality of the system. This behavior was attributed to 
the strong dispersion around the gapless points in the 
quasiparticle excitation spectrum. In order to obtain an- 
alytical results, the authors of Ref. 5 assumed a rel- 
atively small number n c of Landau levels occupied by 
the electrons participating in the superconducting pair- 
ing (n c = Ep/htOc, where Ep is the Fermi energy and 
u> c = eH/m* is the cyclotron frequency). However, this 
assumption is not expected to be quantitatively valid in 
the typical range of fields used in experimentsS^SiSiiiSSiS 7 .. 
On the contrary, the number of occupied Landau lev- 
els n c in the mixed state H c i << H < H C 2 is often 
quite large, typically n c ~ 30 — 270 for borocarbide and 
n c ~ 250 — 4500 for A-15 superconductors. The intention 
of the present work is to numerically evaluate the quasi- 
particle specific heat in the mixed state starting from 
the high-field limit of the Landau level pairing scheme, 
but under more realistic assumptions for the microscopic 
properties of materials studied so that a direct compari- 
son can be made with available experimental data. 

We begin by considering the density of states Af(uj, H) 
in a dirty but homogeneous superconductor in the pres- 
ence of nonmagnetic impurities in high magnetic field. 
In such a superconductor, the coherence length is much 
longer than the effective range of the impurity potential, 
so that under these conditions the order parameter in 
the mixed state A(r) is not substantially affected and still 
forms a perfect vortex lattice. We follow Green's function 
perturbative approach to impurity effects in high mag- 
netic fields developed in our previous workii. Normal 
and anomalous Green's functions for a clean supercon- 
ductor are expanded in terms of the complete set of eigen- 
functions in a magnetic sublattice representation (MSR) . 
In the Landau gauge A = H(—y, 0, 0), the eigenfunction 
^fc 2 ,q,n( r ) belonging to nth. Landau level can be written 
as 
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where £ is the spatial coordinate and k z is the momen- 
tum along the field direction, a = (a, 0) and b = (b x , b y ) 
ar c the unit vectors of the triangular vortex lattice, I = 
y^h/eH is the magnetic length, and V is the volume of 
the system. H n (x) is the Hermite polynomial of order n. 
Quasimomentum q, perpendicular to the direction of the 
magnetic field, is restricted to the first magnetic brillouin 
zone (MBZ) spanned by vectors Qi = {b y /l 2 ,-b x /l 2 ) 
and Q2 = (0, 2a/ 1 2 ) £ In this representation, the "Fourier 
transforms" of superconducting Green's functions in this 
quasimomentum space expressed in the Nambu formal- 
ism can be written as 
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is the quasiparticle excitation spectrum of the super- 
conductor in high magnetic field near the points k z — 

' " This spectrum is 



±k Fn = J 2m(n - hw c (n + 1/2)) /h 2 

calculated within the diagonal approximation 5 ^, where 
only the electrons belonging to mutually degenerate Lan- 
dau Levels at the Fermi surface are involved in the su- 
perconducting pairing. Contributions to the pairing from 
the Landau levels separated by huj c or more are included 
in the renormalization of the effective BCS coupling con- 
stant [g — > g(H, T)]i28i For quasiparticles near the Fermi 
surface (k z ~ fcf n ) it is enough to consider only the 
E n ,p=o bands . The gap function, A nm (q), in the MSR 
can be written as 
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The function A nm (q) turns to zero on the Fermi sur- 
face at the set of points in the MBZ with a strong linear 
dispersion in q. The excitations from other, p / in 
Eq. (J3J, bands are gapped by at least the cyclotron en- 
ergy hu> c and their contribution to a superconductor's 
thermodynamics can be neglected at low temperatures 
(T <C A(T,H) < hui c ). Once the off-diagonal contribu- 
tion is included in the superconducting pairing, the ex- 
citation spectrum cannot be written in the simple form, 
Eq. J2J), and a closed analytic expression for the super- 
conducting Green's function cannot be found. A detailed 
study of the effects of off-diagonal terms on the super- 
conducting state has been pioneered by Norman, Mac- 
Donald and collaborators (see Ref. 30 and references 



therein). Still, when these off-diagonal terms are treated 
analytically within the perturbation theory of Ref. 6, 
the qualitative behavior of the quasiparticle excitations 
at the Fermi surface, as characterized by nodes in the 
MBZ, remains the same. This statement is correct in all 
orders of perturbation theory and therefore is exact as 
long as the perturbative expansion itself is well defined, 
i.e. as long as the magnetic field is larger than some 
critical field H*(T). The critical field H* at T ~ can 
be estimated from the dHvA experiments to be ~ 0.5H C 2 
for A- 15 and ~ Q.2H C 2 for borocarbide superconductors 3 . 
Recent measurements of thermal transport in borocar- 
bides in the mixed stated suggest a strong anisotropy 
of the s-wave order parameter, so that the value for H* 
in these systems can be even lower than the estimate 
obtained from dHvA measurements. Once the magnetic 
field is lowered all the way to H* , the contribution of 
the off-diagonal pairing terms becomes essential 3 ^ and 
gaps start opening up at the Fermi surface signaling the 
crossover to the low-field regime of quasiparticle states 
localized in the cores of widely separated vortices 8 ' 9,30 . 

In the presence of disorder the bare Green's function 
in Eq. @ is dressed via scattering through the diagonal 
(normal) self-energy T, N (iui) and off-diagonal (anoma- 
lous) self-energy E^ n (q, iw)*^ A dressed Green's func- 
tion is obtained by replacing uj with Co and A„„(q) with 
A nn (q) in © where 
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iui = iu — £ (iui) 
A nn (q) = A„„(q) + S^ n (q, iu). 



(5) 



We follow a T-matrix approximation originally devel- 
oped for heavy fermion superconductors^! and adapted 
by us to treat self-consistently impurity scattering at high 
magnetic fielcUi. Within this approximation both weak- 
scattering and strong scattering limits can be treated on 
equal footing. However, we anticipate that the exper- 
imentally determined disorder parameters will put our 
calculation into the weak scattering limit of this the- 
ory with a dilute concentration of impurities. Within 
a T-matrix approximation the self-energies T, N (iui) and 
£^„(q, iu) of the superconducting system are closely re- 
lated to the diagonal (with respect to the magnetic trans- 
lation group basis) T-matrix elements in a single-site ap- 
proximation asii 



sl(q; iw) = 



(iu) = n, < T r "(fc z ,q,w) > Ri 
n i < r^(fe,q,«w) > Ri 



(6) 



where < ... >r,j denotes the average over the im- 
purity positions and n-i is the impurity concentration. 
T^ n (k z , q, ioj) are the coefficients in the T-matrix expan- 
sion over the complete set of MSR eigenstates. The 2x2 
matrix T(r, r'; iui) obeys the Lippmann-Schwinger equa- 
tions 



T(r, r'; iui) = U(t)5(y - r')a z + J dr x U(r) 

xd- z g(r,r 1 ;iw)f(r 1 ,r';iu) (7) 
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where ^-matrix elements are given by Eq. J3J and 
^( r ) = J2i UoS(r — Ri) represents a short-range impu- 
rity potential with the location of impurities Ri taken to 
be randomly distributed everywhere in the sample. The 
scalar scattering amplitude Uq is assumed to be isotropic. 
In high magnetic field we can assume that the scattering 
potential is weak compared to the separation between 
Landau levels, given by hu c . Under these circumstances 
the electrons scatter into the states belonging to the same 
Landau level so that the off-diagonal (with respect to 
Landau level index n) T-matrix elements in Eq. 10 can 
be neglected. Note that this approximation might be 
valid even under more general circumstances, when the 
effective scattering is larger than Tuj c , as is the case in the 
Quantum Hall Effect (QHE) problems 3 ^. However, our 
case, unlike the QHE, also contains strong pairing poten- 
tial and thus the interference between the two is generally 
a rather formidable problem. Following the formalism 
outlined in Ref. 11, Lippmann-Schwinger equations JJJ 
are averaged over the impurity position and from there 



can be solved as 
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where G^(fc z , k, ioS) is a matrix element of a Nambu ma- 
trix (J2J), /„„(k) = A„„(k)/A and V is the volume of the 
system. goes over the entire MBZ while ^2 m is over 
all the occupied Landau levels. Replacing uj with Q and 
A nn (q) with A„ n (q) in (J2J with the help of definitions 
(J5J and |JBJ, Eqs. iJSJ and © can be brought to the form 



w , r Zn=o (rn */4ir 3 k Fn N(0)) J dq(l - V2\f m (q)\ 2 ) u/^U* + |/„„(q)| s 
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where £ = T/A and u = Q/A. N(0) is the normal state 
density of states at the Fermi level. Disorder is character- 
ized with two parameters: T = rii/N(0)Tt = (m/n)EF, 
which measures the concentration of impurities rela- 
tive to the electron density n, and c = 1/ttN(0)Uo, which 
measures the strength of the scattering potential. The 
weak-scattering limit in l|lUfl is approached when c 2 is 
much larger than the second term in the denominator of 
jXQ| , typically when c ~ 1. The strong (i.e. unitary) 
scattering limit is approached when c = 0. The nor- 
mal state inverse scattering rate To is found by taking 
/nn(q) = in (|10f> and letting u> — ► 0. This procedure 
yields 



of disorder is defined as 



r n = 



l + c 2 



(11) 



and establishes the connection between the experimen- 
tally determined disorder parameter To and the param- 
eters r and c in our theory. Equation (|ll)fl is an implicit 
equation from which u = f(u)/A) has to be calculated 
numerically. In fact, once the analytic continuation to 
real frequencies iuj = lu + i5 is performed, this equation 
transforms in a nonlinear system of equations. Finally, 
once u is known the superconducting density of states 
and other thermodynamic quantities can be calculated. 

The superconducting density of states in the presence 
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where Gn(k z , q, iu>) is a Nambu matrix J2J in which the 
replacement lo with Q and A nn (q) with A„„(q) has been 
implemented. Once the analytic continuation to real fre- 
quencies is performed in Eq. I|10[) and Eq. I|12l) and with 
the help of the definition u = Q/A, the density of states 
(|12f) can be obtained as 



N{us,H) _ 1 
N(0) ~ N(0) 



^ 47T 3 fc F „ 
n— 



dq 



V\fnn(q)\ 2 -^ 



(13) 

In deriving this equation, we assume that the impurity 
scattering does not change the q-dependence of the quasi- 
particle excitation spectrum, i.e. A(q) = A/„„(q) is 
given by Eq. (@J. This assumption is shown to be valid 
for the high-field quasiparticle excitations close to the 
gapless points at the Fermi surface while it is less reliable 
for excitations gapped by large A^» We are primarily in- 
terested in the behavior of the specific heat when T — > 
which is governed by the excitations around nodes while 
the contribution of the gapped regions in the MBZ is ex- 
ponentially small. Therefore, this approximation seems a 
modest sacrifice in the quantitative accuracy when faced 
with the overwhelming numerical difficulty in determin- 



5 



ing fully self-consistent A(q) in the presence of the dis- 
order. 

Once the superconducting density of states is obtained 
from Eq. (| 1 3|) the quasiparticle contribution to the spe- 
cific heat in the mixed state at low temperatures can be 
computed asS 

C(T, H) = -J duiNiu;, H) ^ cosh" 2 — (14) 

and then ls (H) = C{H)/T when T -> can be de- 
termined. The most challenging part of our density of 
state and specific heat calculations for realistic super- 
conducting materials in the range of the fields used in 
experimenta 23 i 24 i 25 i 26 i 27 is solving Eq. ijTUIl when the 
number of Landau levels involved in the superconduct- 
ing pairing n c = Ep/huj c is very large. For the A-15 su- 
perconductor VaSi, we estimate n c ~ 241 at H C 2 = 18.5 
Tesla and n c ~ 4470 at H = 1 Tesla (we used an effective 
mass of m* = 1.7m e and Fermi velocity vf = 2.8 x 10 7 
cm/s from Ref. 34 in this estimate). For the borocarbide 
superconductor YNI2B2C, we get n c ~ 25 at H C 2 = 8 
Tesla (from Ref. 24) and n c ~ 200 at H = 1 Tesla, 
obtained using m* — 0.35m e , vf — 2.5 x 10 7 cm/s, a 
normal scattering rate To = 0.53 meV from the dHvA 
experiment of Terashima et al. 3 (these values reproduce 
the mean free path I — 1500 A reported for a clean sam- 
ple by Nohara et al. 25 ). Faced with the overwhelming 
computational difficulty of solving a nonlinear system of 
equations that Eq. 1)10(1 becomes for n c 3> 100 , we limit 
our study to the borocarbide superconductor YM2B2C 
only. 

Figure 1 represents the quasiparticle density of states 
J\f(u>) in the mixed state of the borocarbide superconduc- 
tor YM2B2C computed from Eq. I|13|) and rescaled by 
the normal state density of state N(Q) as a function of 
the reduced energy lu/A. Disorder parameter c in Eq. 
((10JI is chosen to be c = 0.65 which, at the same time, 
determines a value for the second disorder parameter £ 
to be £ = 0.33 if the experimentally determined normal 
state scattering rate of Tq/ A from Ref. 3 is to be repro- 
duced using Eq. I|ll|) . It can be seen in Fig. ^that the 
density of states in the mixed state Af(u}) at low energies 
diminishes as the magnetic field is lowered from H < H C 2 
to H ~ 0.2H c2 . This is a consequence of the depletion 
of gapless or near gapless quasiparticle excitations at the 
Fermi surface at lower magnetic fields. 

The quasiparticle specific heat at low temperatures as 
computed from Eq. (|14fl is, to leading order, a linear 
function of temperature T due to the creation of a finite 
density of states at the Fermi level in Fig. 1. In Fig. 2, 
we plot C(H,T)/T as T -> (i.e. the coefficient of a 
linear-T term in the quasiparticle specific heat) normal- 
ized by the Sommerfeld constant 7 s (i?)/7Ar for YNI2B2C 
as a function of the reduced magnetic field H/H C 2. Full 
circles represent the experimental data of Nohara et al. 
(Ref. 24). We plot two theoretical curves calculated from 
Eq. with c = 0.60 (full line) and c = 0.65 (broken 
line). The second disorder parameter £ is calculated from 
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FIG. 1: Quasiparticle density of states Af(u>) in the mixed 
state rescaled by the normal state density of state N(0) as a 
function of the reduced energy ui /A computed from Eqs. (11 0> , 
1111 and 1131 1 for the borocarbide superconductor YM2B2C 
when disorder parameter c = 0.65. We have used experimen- 
tally determined values H C 2 = 8 Tesla and T c — 15.4 Kelvin 
from Ref. 24 as well as m* = 0.35m e , To = 0.53 meV and 
A = 2.3 meV from the dHvA measurement of Terashima et. 
al. (Ref. 3). 




FIG. 2: The coefficient of the linear-T term in the quasipar- 
ticle specific heat at low temperature normalized by the Som- 
merfeld constant r y(H)/"/N for the borocarbide superconduc- 
tor YNi2B2C as a function of reduced magnetic field H/H C 2- 
Full circles represent experimental data of Nohara et al. (Ref. 
24), while lines represent our theoretical curves calculated 
from Eq. O with the help of Eqs. JTOJ, and (EHJ for 
two values of the disorder parameter c (see text). 



Eq. (|11|) where the experimentally determined value for 
the normal state scattering rate To = 0.53 meV from 
Ref. 3 is used. The values of the other physical quanti- 
ties needed in our theory, effective mass m* = 0.35, BCS 
gap A = 2.3 me V and upp er critical field H C 2 = 8 Tesla 
in A(H) — Ay/1 — H/H C 2, are taken from experiments 
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of Nohara et al. (Ref. 24) and/or Terashima et. al (Ref. 
3). Fig. 2 demonstrates that the specific heat coefficient 
computed from the theory presented in this paper ex- 
hibits a nonlinear dependence on the magnetic field H, 
i.e., j s (H) ~ H 37 for c = 0.60 and j s (H) ~ H 0A6 for 
c = 0.65. This is in contrast to the linear H dependence 
predicted for a fully gapped s-wave superconductor in 
the mixed state as T — ► 0. We attribute this nonlinear 
behavior to the creation of coherent, gapless quasiparti- 
cle excitations at the Fermi surface in the mixed state 
of an extreme type-II superconductor at high magnetic 
field H such that H* < H < iJ c2 . The estimated critical 
field H* for pure YNi 2 B 2 C is - (0.15 - 0.2)H c2 & How- 
ever, this estimate depends on the value of the s-wave 
gap function and can be much smaller than 0.15iJ C 2 if 
the value of the minimum gap in the strongly anisotropic 
s-wave case is significantly different from the accepted 
BCS valueit Furthermore, it seems that the unusual 



behavior of the specific heat coefficient J S (H) is also a 
consequence of disorder present in the superconducting 
system. The finite amount of impurities in the system 
leads to the creation of a finite density of states at the 
Fermi level Af(0), in contrast to the perfectly clean su- 
perconductor where Miui) exhibits an algebraic behavior 
at low energiesii. 

In summary, we have computed the quasiparticle con- 
tribution to the specific heat C(H, T) within the T- 
matrix approximation for the self-energies and found a 
nonlinear behavior of j s (H) = C(H,T)/T when T —> 0. 
In the range of magnetic fields where our theory is ap- 
plicable H* < H < H c2 , the calculated 7«(-ff) closely 
resembles the experimental data for the borocarbide su- 
perconductor YNi 2 B 2 Ci2£ 
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